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According to a continuous medium theory, in very thin ferroelectric films with real metallic 
electrodes (or dead layers near the electrodes) the domain structure reduces to sinusoidal distribution 
of ferroelectric polarization. Such a sinusoidal structure was considered in 1980s for para-ferroelectric 
phase transition in a capacitor with dead layers near electrode. We give a review of this theory and 
its further development for the case of real metallic electrodes. The goal of the general theory is to 
consistently interpret the experimental data in very thin films with real metallic electrodes. This is 
illustrated on a recent experimental data for 5-30 nm BaTi03 films with SrRuOs/SrTiOs electrodes. 
The screening length by real metallic electrodes is very small (< lA), but it has a profound effect 
on ferroelectric properties and its phase behavior. This general theory also allows to formulate the 
important open problems and show paths towards their solution. In particular, this is a problem of 
finding parameters of the system, which can sustain the ferroelectric memory over a desired lifetime. 

PACS numbers: 77.80.Dj, 77.55. +f, 77.22.Ej 
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I. INTRODUCTION 

Apparently, the continuous medium theory of domain 
structures in ferroelectrics (FE) becomes even more im- 
portant nowadays than in prior years, when the studies 
concentrated on properties of the bulk samples. One of 
the reasons is that the phenomenological domain the- 
ory, established decades agoi^^, consistently captures 
the electrostatic and strain effects playing crucial role 
in very thin high quality FE films available now, and al- 
lows to gain an indispensable qualitative insight into the 
phenomena at hand. One expects that the theory should 
break down in atomically thin ferroelectrics, where the 
calculations should be done from the first principles. This 
is formally true, yet the theory applies to most systems 
of interest that are not just a few monolayers thin, since 
the atomic length is so small. And the FE systems of 
most interest to e.g. memory applications would not be 
that thin. In fact, they should be thick enough to sup- 
port a memory of practical interest, and we reach this 
conclusion safely within the domain of applicability of 
the continuum theory, as described below. The present 
review has a restricted goal: it is aimed at explaining 
what the continuous medium theory tells us about the 
properties of thin films with real metallic electrodes, es- 
pecially as far as domains and ferroelectric memory are 
concerned. We shall discuss not only the theory but also 
some recent data, and we hope that as a result of this 
discussion, the important role of the phenomenological 
theory as well as necessary further developments of this 
theory, will become rather evident 

The principal point of the review is that some old the- 
oretical results about specific features of domain struc- 
tures near a second order FE transition, which were 
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mainly of academic interest, gained practical importance 
now. Let us recall that close to a second order tran- 
sition the ferroelectric domain structure cannot remain 
the same as it were far from the transition, where do- 
main walls thickness is smaller than the width of the 
domains. The reason it that the domain wall thickness 
tends to infinity when the temperature approaches the 
critical temperature T c , T — > T c , while the period of the 
standard domain structure remains almost temperature 
independent. As a result, some distance away from T c 
the width of the domain walls becomes comparable with 
the width of the domains, i.e. the spatial distribution of 
the polarization, which represents the domain structure, 
becomes pretty smooth ("sinusoidal"). One can picture 
this change as elimination of higher Fourier harmonics 
of the initially staggered distribution of the polarization 
across the domain structure. The main harmonic is, ob- 
viously, the last to disappear since it defines the main 
sluggish feature of the domain structure: its period. Dis- 
appearance of the main harmonic is the phase transi- 
tion into the paraelectric phase, which occurs somewhere 
at temperature, T < T c . For macroscopic samples with 
thickness I 3> d at , where d a t is the characteristic inter- 
atomic distance, the elimination of the higher harmonics 
is in fact quite sharp: all the above described transfor- 
mation occurs within a tiny, experimentally inaccessible 
temperature interval near T — T c . The thinner the film 
the broader this interval, and when the film thickness 
approaches the unit cell size, it may become quite im- 
portant. Indeed, the period of the domain structure is 
about y/daJ, Ref4 (see also^i^ and our discussion be- 
low), while far from T c the domain wall thickness can 
be roughly estim ated as d at for the order-disorder tran- 
sitions and d a t\/T a t/T c for the displacive ones, where 
T at ~ 10 4 — 10 5 K is the so called characteristic "atomic" 
temperature. These estimates are fairly rough, for exam- 
ple in BaTiOs, which is usually considered a displacive 
ferroelectric, the 180— degree domain walls parallel to the 
direction (100) and the equivalent ones are actually very 
narrow with the width of about d a t at room temperature. 
However, the estimates clearly indicate that in very thin 
films, with thicknesses approaching d at and for displacive 
systems, and maybe even in thicker films not far from 
the phase transition, the domain width and the domain 
wall thickness become comparable, so that the domain 
structure is expected to be "sinusoidal" . 

It is worth mentioning that the basic physics of an 
effect of depolarizing field and domains on the phase 
transition was well understood since the pioneering work 
of Kanzig's group 9 in 1950s, but the case of the "sinu- 
soidal" domain walls was first considered by Chensky and 
Tarasenko (ChT) some 30 years later—. Historical as- 
pects of the topic are described in more detail in Ap- 
pendix A. ChT were interested in the second order phase 
transitions in ferroelectric films in both non-electroded 
and electroded samples with a dielectric ("dead") layer 
separating the electrodes and the film. In the FE ca- 
pacitors with metallic electrodes, the role of the "dead 



layers" is played by the metallic electrode interfacial re- 
gions over the Thomas-Fermi screening length. Within 
the continuous medium theory the mathematical analogy 
between the two cases is practically exact (see below). 
Therefore, the model considered by ChT was, in fact, 
very general. The results attracted little attention so far, 
since the sinusoidal domains were expected to exist in 
a minute temperature interval near phase transition in 
the systems experimentally studied at the time. Now, 
however, the situation is quite different. In ferroelectric 
films of atomic thickness, the sinusoidal regime has been 
experimentally found to exist in a temperature interval 
of about 100K near the transition 1 ^. The approach it- 
self, new results, and the relevance of the theory to the 
present experimental studies of nanoferroelectrics need to 
be exposed to the ferroelectrics community. The theory 
becomes of practical interest, and new problems need to 
be addressed. We discuss, in particular, the problem of 
a (meta)stable ferroelectric memory. While most people 
are currently concentrating on the problem of a critical 
thickness for ferroelectricity, which is shown to approach 
the atomic limit, the real practical problem is that of 
a ferroelectric memory, i.e. a (meta)stable state with a 
net ferroelectric polarization in the film. These are two 
different problems, and we are trying to highlight this 
important distinction. 

Below, in Sec. II we present a simple overview of the 
ChT treatment of paraelectric-ferroelectric transition in 
a film with the ideal metallic electrodes and the dead 
layers. At the end of this Section, we demonstrate sim- 
ilarity of the cases of real metallic electrodes and the 
dead layer treated in 1 ^. In very thin strained (and, as 
a consequence, uniaxial) films considered here, where a 
domain structure is close to the sinusoidal one, we show 
the distribution of the polarization, which is reminiscent 
of the one for the domains with staggered distribution of 
polarization (the so-called Kittel structure^.) 

In Sees. Ill, IV we consider the properties of sinusoidal 
domain structures, mainly their response to an external 
field. The main results of these Sections present a fur- 
ther development of ChT theory. The algebra involved 
is relatively simple in cases when the Landau-Ginzburg- 
Devonshire (LGD) free energy contains the powers of 
polarization P up to P 4 terms only (in the bulk or in 
the case of metallic electrodes this would correspond to 
a second order transition far from a tricritical point), 
or just P 2 and P 6 terms (describing a tricritical tran- 
sition in the same conditions), while in a more general 
case of a second order and/or the first order transitions 
it is fairly involved. We limit ourselves here by the two 
above cases, which illustrate the main results well enough 
to enable us to discuss the available experimental data 
in Sec.V. Importantly, irrespective of the conductivity 
type of electrodes, metallic or semiconducting, the im- 
perfect screening of bound charge of the FE films leads 
to their splitting into domains 12 . The results— are used 
to discuss the experiments by the Noh's team on BaTiOs 
films with SrRuO electrodes on SrTi03 substrat o 14 ' 15 i 16 
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in Sec.V. There, we conclude that their 5nm thick sample 
is well above the critical thickness for ferroelectricity, but, 
even by very optimistic estimates, it is below what can 
be called the critical thickness for the (meta)stable ferro- 
electric memory with a reasonable lifetime (e.g. in excess 
of 10 3 }i£. Unfortunately, at present the theory cannot 
give the realistic number for this thickness. Throughout 
the review, we omit the so-called additional boundary 
conditions (ABCs), which include parameters of the sur- 
faces/interfaces and are often taken into account together 
with conventional electrostatic boundary condition a 18 ' 19 . 
For the experimental system that we discuss here, these 
additional boundary conditions appear inessential, and 
more details about the additional boundary conditions 
can be found in Appendix B. 



II. LOSS OF STABILITY OF THE 
PARAELECTRIC PHASE 



U=ax 2 




FIG. 1: Schematic of the potential energy U(x) = ax 2 as a 
function of coordinate x for various a: a > 0, the state with 
x = is stable; a = 0, point of stability loss of the state with 
x = 0; a < 0, stability at x = is lost. 



The transition point is simultaneously the point of a 
stability loss of the two phases for any second order phase 
transition. We shall consider second order transitions 
from paraelectric phase in an electroded ferroelectric film. 
This transition can be either into a single domain or a 
multidomain state. To decide what choice is being real- 
ized (and at what conditions) the best way is to consider 
the loss of stability of the phase whose state we know: 
the paraelectric phase. 



dead/screening layer 



FE 



dead/screening layer 



A. The method 



Let us explain how to determine when the loss of sta- 
bility occurs by using an elementary example. We show 
three forms of a potential energy U (x) = ax 2 versus 
some coordinate x in Fig. [TJ Obviously, the state with 
x = is stable for a > 0, for a < this state is un- 
stable, and a — (the horizontal line) corresponds to 
the point where the stability is lost. The condition of 
equilibrium dU/dx = has only the trivial solutions for 
a ^ 0, but an infinite number of nontrivial solutions for 
a = 0. For the potential energy of a more general form, 
U (x) = ax 2 + bx 4 (b > 0), we have the same point of the 
stability loss, but dU/dx = has the trivial solution only 
for a > 0. However, the linearized form of U near the ori- 
gin x — is the same as before. We have, therefore, the 
method of finding the conditions of instability of a given 
state (or phase) : one has to linearize the equations whose 
solution represents the state in question and look for the 
conditions when this system of equations has a nontrivial 
solution. 

We shall consider a uniaxial ferroelectric with the fer- 
roelectric axis perpendicular to the film plane, see Fig. [5] 
This is the case, for example, of BaTiC>3 (BTO) films 
compressively strained on SrRu0 3 /SrTi0 3 (SRO/STO) 
substrate. The strained cubic crystal behaves as a uni- 
axial ferroelectric with second- or weak first-order phase 



FIG. 2: Schematic of the ferroelectric thin film with either 
dead layers with thickness d/2 separating it from the elec- 
trodes, or the real electrodes with Thomas-Fermi screening 
length A. 

transition. We restrict ourselves here to the case of a 
second order ferroelectric transition in the bulk, i.e. the 
Landau-Ginzburg-Devonshire (LGD) free energy density 
is assumed to have the form: 

F LGD (P) = ^P 2 + ^ + ... + ^D t3 (V ±t P z )(V ±J P z ) 
+ l - 71 { dz P z ) 2 + l - K p 2 b (1) 

where P z is the ferroelectric (switchable) component of 
polarization along the polar axis z, P z t is the nonfer- 
roelectric part of the polarization along the z-axis with 
K directly related to the so-called background dielec- 
tric constant (see the Methodological Note below), A = 
A' (T - T c ) and A', B, D tf , 77, k - const, V x = (d x ,d y ) 
the gradient operator in the plane (x — y) of the film. 
To simplify the formulas, we shall assume that the gra- 
dient coefficient is L\j = DSij, where dij is the Kro- 
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necker symbol. Here, the Landau coefficients are renor- 
malized by the strain produced by the lattice misfit with 
the substrate Unr^, see the Appendix D (for BTO on 
STO u m = -2.2%). 

Importantly, we do not include the striction terms that 
couple strain to P 2 terms in ([I]) because they renormalizc 
the BP^ term and this does not affect the loss of stabil- 
ity of symmetric phase in a system like BTO on STO 
(see below). The point is that, as we said above, finding 
the loss of stability boils down to solving a linear prob- 
lem, while this term is nonlinear and does not affect the 
problem. The statements to the contrary by Pertsev and 
Kohlstcdt in Refill, that the elastic coupling defines the 
point of the stability loss and the critical thickness of the 
FE films, are qualitatively incorrect, as pointed out in 
Refi22. The above form of Flqd yields the equation of 
state for the ferroelectric via standard relation: 



E 



5F, 



LGD 



SP Z 



AP Z + BP* - DV\P Z - nd 2 z P z , (2) 



E fz = kP : _ 



zb • 



(3) 



which means that the local electric field is in one-to-one 
correspondence with the local ferro- and nonferroelectric 
components of the polarization. 



B. Loss of stability with respect to homogeneous 
polarization 

Since we do not know what type of solutions corre- 
sponds to the horizontal line in Fig. [TJ we should check 
all the possibilities. We start by trying the simplest one, 
a homogeneous change of polarization. This situation is 
described by the linearized equation of state for a homo- 
geneous polarization P z h (and the nonferroelectric part 
KP zb ) : 



E 



fz 



AP 7 . 



kP* 



(4) 



where Ef z is the homogeneous field in the ferroelectric. 
The dead layer is a linear dielectric, which is described 
by an equation analogous to hence it is not necessary 
to consider explicitly the polarization of the dead layer, 
and we shall describe it by its dielectric constant e e . The 
boundary condition at the interface of the ferroelectric 
(continuity of the normal component of the displacement 
vector) then reads 



E fz + An (P zh + P zh ) = e e E dz 



(5) 



where Ed z is the field in the dead layer. We can rewrite 
the above equation as 



£bEf z + AnP zh = e e E dz , 



(6) 



where e b is the "background dielectric constant", e b — 
1 + An/ n. There is also a condition of a zero applied bias 
(short circuiting, {7 = 0): 



E fz l + E dz d = 0. 



(7) 



We find from the above the depolarizing field in the FE 
film: 



Ef z - 



E, 



And 



-,P: 



AndP 7 



ej + e b d 
4tt/ 

-Pzh 



zh 



6,7 



4ttR 



e e l + e b d 



(8) 
(9) 



and we see that the field in the dead layer is much 
stronger than the field in the FE film, \E dz \ ^> \Ef z \ be- 
cause is does not contain the small parameter d/l <C 1, 
and it is realistic to assume e b ~ e e . Importantly, the field 
is concentrated in the dead layers, yet since the layers are 
very thin, its contribution to the total free energy appears 
to be negligible (see below). The system of homogeneous 
equations (J4j) , ((8j> for Ef z and Pq has nontrivial solutions 
only when 



.4 



And 



e e l + e b d 



And 



(10) 



This would define the temperature of the stability loss of 
the paraelectric phase with regards to a homogeneous po- 
larization. It would happen at a lower temperature than 
the loss of stability of the paraelectric phase in a very 
thick slab (I — > oo) or in absence of a dead layer and ideal 
metallic electrodes (A = 0, T = T c ). Recall that absence 
of nontrivial solutions at values of A which do not sat- 
isfy Eq. (JTOj) means cither stability of the solution P = 
(the paraelectric phase) or instability of this solution (see 
Fig.l). and, therefore, stability of another homogeneous 
solution with P / 0, which has to be found from the 
full, nonlinearized, equation of state. Let us emphasize 
that this is not necessarily the absolutely stable state 
but might be a metastable state because we did not con- 
sider competing inhomogeneous states which might have 
a lower energy. Also the temperature defined by Eq. (|10|) 
is not, in general, the temperature at which the paraelec- 
tric phase ceases to exist: one has to consider a loss of 
stability with respect to an inhomogeneous polarization 
as well, in order to determine which one sets in first. 

Now, consider the case of a finite bias voltage U. In 
this case, the equilibrium polarization in the paraelectric 
phase is non-zero, and so are the electric fields in the 
ferroelectric and the dead layer. One finds 

lE fz + dE dz = U 

and, instead of Eq. (JH), we obtain the field in the FE film: 

And „ 



Efz — Eq 
Eq 



e e l + e b d 

U 

I + e b d/e e I 



(11) 
(12) 



where Eq is the external field. The expression for Ef z 
(fTTj) is usually said to mean that the field in the FE film 
is equal to a sum of the external Eq and the depolarizing 
fields. Substituting Eq. (TTTj) into Eq. ©, we obtain 



A 



And 



e e l + e b d 



P Q + BP3 = E , (13) 
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where Po stands for the equilibrium homogeneous polar- 
ization. 

Let us check if this equilibrium state can lose its sta- 
bility with respect to a homogeneous fluctuation of po- 
larization. In textbooks on ferroelectricity, it is shown 
that when d = (an ideal metallic electrodes) there is 
no loss of stability, and the second order transition gets 
smeared out by the external field. Lets check that the ho- 
mogeneous transition would also be smeared out in the 
case of the dead layers. Checking for the stability loss 
with respect to a homogeneous polarization, we linearize 
the equation of state ^ near P = P , P/ z o, where Ef z0 
is the electric field in the ferroelectric in the state with 
P = Po, and obtain 



(A + 3PP 2 ) P' = E'j 



(14) 



Po and Ei = Ef z —Ef Z Q. We obtain the 



where P' = P 

same equations for the fluctuation of the polarization P' 
and the field Ej z as for the short-circuited case with the 

only difference that A should be replaced by A = A + 
3PP 2 . We have already found [cf. Eq. $TU})} that the loss 
of stability with respect to a homogeneous polarization 
should take place if 



A = A + 3-BPq 



And 



e e l + e b d 



(15) 



This is impossible, however, since from Eq. |T 



A 



And 



e e l + e b d 



A 



And 9 

— + 3PP 2 

ej + e b d 



dP Q 
dE 



(16) 

is always positive. Indeed, dP^/dE^ can be considered 
as the differential susceptibility of a ferroelectric in an 
external electric field and with the ideal metallic elec- 
trodes, but with the coefficient before the P 2 term equal 
to A + Aird/ (e e l + ebd). We have already mentioned that 
in the case of the ideal metallic electrodes there is no 
loss of stability and the phase transition is smeared, i.e. 
dPo/dEo is always positive under a finite bias voltage. 
We have proven that in a ferroelectric with the dead lay- 
ers (or real metallic electrodes) the phase transition with 
respect to the homogeneous polarization is smeared out 
by an external electric field. As was mentioned above, 
this does not mean that there is no phase transition. In- 
deed, there is a possibility that in an external field the 
film may be losing stability with respect to inhomoge- 
neous polarization. If this takes place, a second order 
(unsmeared) phase transition occurs This is why we shall 
study a loss of stability with respect to the inhomoge- 
neous polarization in a general case of applied external 
bias voltage. 



C. Loss of stability with respect to inhomogeneous 
polarization 

Following the same general procedure, we should be 
looking for a first occurrence of a nontrivial solution to 



a linearized equation of state that allows for inhomoge- 
neous polarization: 



E z =AP z -DV\P z -i 1 d 2 z P z . 



(17) 



Since an inhomogeneous polarization P z produces an 
electric field both across the film (z— axis) and in perpen- 
dicular directions, we should add an equation of state for 
the two other polarization components. We have: 



- 1 



± = 



An 



-E\ 



(18) 



where E± is the in-plane electric field, e± the in-plane 
dielectric constant. The solution of (fT7|) should also obey 
the equations of electrostatics: 



Vx£ = 0, 



V- (e b E + AnP) = 0. 



(19) 



(20) 



Now, we should get inventive. As far as a reasonable solu- 
tion of equations (fTT|) , (fl"9]) , (|20|) is concerned, they have a 
vast set (continuous infinity) of types of solutions and we 
have to single out only one type of solutions from them 
that gets realized. We shall apply a physical common 
sense to restrict the types of the solutions considered. 
First, the prospective solutions should be periodic in the 
film plane, since we deal with a large area (practically 
infinite) film. Any periodic function can be presented as 
a sum of the sinusoidal functions. Consider the structure 
that is periodic in the x-direction. Then, 



P z (x,z) = ^Pn (z)cosk n x, 



(21) 



where k n — nki, P n [z) are the unknown functions and 
ki is the unknown wavenumber. We have to perform the 
stability check for arbitrary k±. The solution given by 
Eq. (f2Tj) is very complicated, and in reality one does not 
need to keep all harmonics. It is sufficient to consider the 
critical harmonic with a one-dimensional periodicity 



Pz — o.f (z) cos kx, 



(22) 



where a and k are the unknown numbers, / (z) the un- 
known function. It adds nothing new if we consider, 
e.g., cos k\x cos k^y, because the products of trigonomet- 
ric functions are the linear combinations of trigonometric 
functions of other (sum and difference) arguments. Look- 
ing among the forms given by Eq. (|22[) . we have already 
tremendously simplified the problem: for every k we will 
have a system of ordinary differential equations, instead 
of a system of partial differential equations. Further, we 
can identify / (z) too. First of all, this function should 
be as smooth as possible. Indeed, the bound charge as- 
sociated with the inhomogeneous ferroelectric polariza- 
tion p b = — V • P = a(df/dz) coskx should be as small 
as possible to reduce the electric field energy. It is un- 
reasonable to have / (z) = const, the bound charge at 
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the ferroelectric-dead layer interface will be proportional 
to this constant. More reasonable would be to make it 
nearly constant deep inside the ferroelectric yet smaller 
at the interface. Taking into account that it should be a 
solution of a differential equation with the constant co- 
efficients, we are left with an ansatz: 



P z = a cos kx cos qz, 



(23) 



where q is another unknown parameter, which we ex- 
pect to be much smaller than k. A critical reader could 
study other options. It may seem formally possible to try 
coshqz, for example, but this is a bad choice that fails 
to satisfy the boundary conditions (see below) , which are 
possible to satisfy with Eq. (|23p. 

Using Eq. (I17[) . one obtains for the ansatz (|23| : 



Efg (x,z) = (A + Dk 2 + r)q 2 ) a cos kx cos qz. (24) 
The electrostatics equation (fT9|) yields 



d z Ef x = d x E 



(25) 



and one finds that 

Ef x (x, z) = (A + Dk 2 + r]q 2 ) a sin kx sin qz. (26) 

Thus, we have found the electric field in the ferroelec- 
tric assuming that we have a "polarization wave" . This 
assumption is consistent if this field satisfies the second 
electrostatic equation (|20| . or 



e b d z E fz + e ± 8 x E fx + 4ird z P z = 0. 



(27) 



We can easily get rid of Ef x in the above equation by 
differentiating it with respect to z and using Eq. (|25p to 
obtain 



(e b d 2 + e ± d 2 ) E fz + A7rd 2 P z = 0. 



(28) 



Substituting here the ansatz (|24p . we find the value of 
A, at which the non-trivial solutions are possible as a 
function of the parameters k and q: 



Anq 2 



(bq z 



e±k 2 



Dk 2 



(29) 



It is still premature to look for a maximum of this func- 
tion (a point, where the instability sets in first), because 
we have also (i) to find the fields from the electrostatic 
equations in the dead layer and (ii) to satisfy the bound- 
ary conditions. Given that the field in the dead layer has 
to have its z component depending on x as cos kx, we 
put Ed z = C (z) cos kx and, therefore, E dx = tp (z) sin kx, 
where ip' (z) — —kC,{z). From V-E = 0, one has 
£' + kip = or %p" — k 2 ip = 0. Since the transversal 
component of the electric field is zero at an ideal metallic 
surface, we find that if> (z) = sinh k (z — 1/2 — d/2) and 

E dx = m sinh k (z-l/2- d/2) sin kx, (30) 
E dz = cosh k (z-l/2 -d/2) cos kx, (31) 



for z > 1/2. The two electrostatic boundary conditions, 
(i) the continuity of the in-plane electric field component 
and (ii) the continuity of the displacement, Eq. ([5]), at 
the ferroelectric-dead layer interface (z = 1/2) read 

a~ (A + Dk 2 + riq 2 ) sin ^ =* s inh^, (32) 



a (A + 47r/e b + Dk 2 



2\ I 1 
m ) cos — 



e e T kd 
— V cosh — . 
£6 2 



(33) 

Excluding A from equations (|32p . (|3"3"f . with the help of 
Eq. (|29|) . one finds that this system has a non-trivial so- 
lutions for a and \& if 



q ql k kd 
— tan — = — tanh — . 
e± 2 e e 2 



(34) 



Since we are interested in the case of thin dead layers, 
we suppose that kd -C 1. Then, Eq. (|3~4"1) simplifies to 



gtan 



ql e\k 2 d 
~2 ~ 2e R 



(35) 



This allows one to exclude k from (f3"5]) and (f2"9"| and re- 
duce, for the case of thin dead layer, i.e. for I ^> deb/e e , 
dr)e±/ (De e ) , to:. 



T 47TC? . . 

A = -jt{y), 
e 2 D 



d 2 ' 



t(y) = V (cot y + r tan y) , (36) 
y = ql/2. (37) 



Now, we can find at what critical thickness of the dead 
layer the first nontrivial solution appears with q =/= 0. 
Consider first the case y = ql/2 <C 1, where the function 
t(y) can be expanded as 



t 



y 



l 

45 



y 



(38) 



There appears a minimum at y ^ for r < 1/3, so that 
the minimal value of \A\ is found at q ^ for the total 
thickness of the two dead layers exceeding the critical 
thickness dm'- 



This solution is: 



d > d rn = 



1 / d- 
l \ d r , 



3D 



1/2 



5n d — d ri 
~T P Dl 



1/2 



(39) 

(40) 
(41) 



We see that for d < d m the function A (q) is monotonous 
and has its maximum value at q = 0, i.e. for very thin 
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dead layers the phase transition occurs into the homoge- 
neously polarized state. In the case d > d m , the phase 
transition occurs into an inhomogeneously polarized state 
with k — k c given by Eq. (|4"Tj) , i.e. the "wave vector" of 
the sinusoidal polarization "wave" increases very rapidly 
from k = at d = d m to a large value and ql/2 rapidly 
becomes on the order of unity. The corresponding tran- 
sition temperature is found from 



\A\ 




5(d— d„ 
4Z~ 



(42) 



i.e. it occurs just slightly above the transition into 
the homogeneous state, where \Ah\ = &wd/e e l. For 
BaTiOa/SrRuOs the estimate gives a small value d m — 
0.3A, hence for BTO/SRO such a regime is only of an 
academic interest, d ~ 5d m there (see Sec. VI). However, 
it seems to be very interesting to look for systems where 
d m may be large, since those systems could perhaps sus- 
tain the homogeneous polarization and, therefore, are 
most interesting for memory applications. 

In the opposite limiting case, d ^> d m , which is a typi- 
cal real situation, ql increases and approaches 7r, so that 
the value of tan & becomes very large. We suppose that 
in this case e±k 2 3> e&g 2 and shall check later that this 
assumption is justified. Then, one can omit ebq 2 in the 
denominator of the first term on r.h.s. of Eq. (|2"9"|) . Also, 
the last term can be omitted, because, as we shall see 
below, the maximum of A corresponds to k 3> q and 
the coefficients D and rj are normally of the same order 
of magnitude. Indeed, the maximum of the simplified 
expression corresponds to 



1/4 



/ 4tt 3 



\e±DP 



1/4 



since, as we argued above, 



q c ~ n/l. 



(43) 



(44) 



In this regime e 1 / 2 k c /q c = (4e ± /7r) 1/4 l^/D 1 / 4 . Given 
that usually D is on the order of magnitude of d 2 t , where 
d a t is the typical interatomic distance (cf. Appendix E) 

and (4ex/7r) 1 ^ 4 is on order of unity, in all real cases one 

has approximately: e 1 j 2 k c lq c ~ (l/dat) 1 ^ 2 , i.e. our as- 
sumption that e±k 2 e^q 2 for d S> d m is justified for 
films thicker than several monolayers. 

We see that with increasing thickness of the dead layer 
d the period of the sinusoidal domain structure rapidly 
decreases and then saturates at a value independent of 
d. Since at d - d m = d rn both Eq. gj) and Eq. ([43]) 
give similar values for k c (difference is about 20%), the 
two simple formulas are sufficient to describe the function 
k c (d) with a reasonable accuracy at all values of d. From 
Eqs. (|i3]>. ([44"f . and (|29j). one sees that at d - d m > d m 
the loss of stability of the paraelectric phase with respect 
to formation of the domain structure takes place at 



The relation (|29|) allows to simplify the expressions 
for the inhomogeneous fields in the ferroelectric film 
(|24| . ([26|) . which we rewrite, using e±k 2 ^> ebq 2 , as: 



Efz (x, z) 



4tt 3 



e± (kl)' 



■ cos kx cos qz, (46) 



4tt 2 



Ef x (x 1 z) — a — — sin kxsinqz. 

6 | rZl 



(47) 



Obviously, both fields are much smaller than the magni- 
tude of the polarization wave: since kl ~ \Jl/d a t S> 1, 



then 



Efx 



< a. 

In comparison, the field in the dead layers in the stud- 
ied case fed <§; 1 is, near FE interface: 



E, 



dz 



i 

x,z = - 



E dx x, z 



a — -7r— cos kx 

e\_k l ld 

47T 1 /2 jD l/2 
11 VST 

e_L d 

4tt 2 . , 
a — — sin kx, 
e±kl 



■ cos kx. 



(48) 
(49) 



hence, E& z is much larger than any other field compo- 
nent in the system. As we can see from (|4"8"|) . there is no 
small parameter like d at /l or d/l in the z— component of 
the field in the dead layer Ed z (the only parameter there 
is d at /d, d at ~ D 1 / 2 [see discussion after Eq.(|33])]. which 
is less than unity, yet much larger than the above men- 
tioned other small parameters in the problem.) This is 
reminiscent of the behavior that we saw above for the case 
of the homogeneous polarization. These results allow us 
to construct the maps of the polarization P z (|23|) and 
P x = (e± — l)Ef x /Air characteristic of the sinusoidal po- 
larization distribution, shown in Fig. [3J The map shows 
the narrow sinusoidal domains with well ordered polar- 
ization distribution. The stray electric field component 
along the film near its surface, E x , results in the pres- 
ence of perpendicular to the easy z— axis component of 
the polarization, P x , in the narrow sinusoidal domains. 
The pattern is similar to the one for the structure with 
abrupt domains (the well-known Kittel structure with 
thin domain walk 2 .) 



D. Metallic electrode 

We would like to compare now the ferroelectric film 
with the dead layers to the one of real metallic electrodes 
with Thomas-Fermi screening, in order to see the analogy 
between the two cases. Obviously, the inhomogeneous 
polarization in the form (|23p will produce the electric 
field and the electrostatic potential oc cos kx, so that the 
electrostatic potential will be cf>{z, x) = &(z) coskx. The 
Poisson equation then takes the form 



.4 



-2Dki 



(45) 



(d 2 - k 2 ) $ = $/A 2 , 



(50) 



8 



/=2nm 

V/////////////A *~. 




x/l 



FIG. 3: The polarization distribution in the BTO/STO film 
with thickness I = 2nm in sinusoidal domain regime. It clearly 
shows that the electric field component along the film near 
its surface results in the presence of perpendicular to the easy 
z— axis component of the polarization in the narrow sinusoidal 
domains. The pattern is similar to the one for the structure 
with abrupt narrow domains (the Kittel structure). 



that has the solution $(z) = — <&(— z) : 

-k(z-1/2) 



$ = pie 
k = \lk 2 + X~ 



z > 1/2, 



(51) 
(52) 



With this solution, the boundary conditions, found anal- 
ogously to ([3"2 j) .(|3"5 ]) . give: 



k nr 
a- (A- 

q 



Dk 2 

a (A + Air/et, + Dk 2 



sin— = kpi, 



(53) 



?\ $ e e 

r/g J cos — = — Kpi. (54) 

4t 6fo 



Excluding A from these equations with the help of 
Eq. (f29l . one arrives at^: 



q , qi 

— tan — = 

e± 2 



; Vfc 2 + A" 



(55) 



Considering the case of interest to us, k\ <C 1, and com- 
paring this with p5|) . we see that for e±k 2 ^> q 2 the 
case of realistic metallic electrodes and the simple dead 
layer model are mathematically equivalent when A = d/2. 
Since d/2 is the thickness of the dead layer, the result 
suggests the numerical equivalence of the dead layer and 
the Thomas-Fermi screening length when e±k 2 3> q 2 ■ 



wave vector equal to k c in addition to a spatially homo- 
geneous part p. There will certainly be a single sinu- 
soidal wave, at least close to the transition, if there is an 
anisotropy in the (x — y) plane. In this review, we discuss 
this latter case only. 



A. Range of existence 

Consider first a zero external field, E = 0, when the 
homogeneous part of the polarization is zero both in 
the paraelectric and the ferroelectric phases, i.e. close 
to the transition there is a single sinusoid and p = 0. 
Some distance away from the transition, the single si- 
nusoidal approximation (|23[) breaks down. Indeed, be- 
cause of the cubic term in Eq. ((2]), the polarization of 
the form P — a cos kx cos qz suggests a presence of an 
electric field, which depends on the coordinates not only 
as cos kx cos qz, but also as cos3fcxcos<7z, cos/ca;cos3gz, 
cos 3kx cos 3qz. Such a field induces polarization contain- 
ing even higher spatial harmonics, and so forth ad infini- 
tum. Let us estimate when these higher contributions 
to the field and the higher harmonics in the polarization 
may be neglected. Suppose that the main part of the po- 
larization is the simplest polarization wave (|2"3"]) . When 
the polarization is known, one can calculate the electric 
field produced by this polarization without referring to 
the equation of state. In our case, it is convenient to use 
the relation between the polarization and the field at the 
transition point given by Eq. (|28|) . which yields: 



E 



fz 



Anq 2 

= — t<2 cos k c x cos q c z 

c±k 2 c 

= — Dk 2 a cos k c x cos q c z, 



(56) 



since e±k 2 q 2 . It is impossible to satisfy exactly 
the equation of state ([2]) with such a field and polar- 
ization, and we have to consider the corrections to the 
field E z (x, z): E z — Ef z + E z . Let us try P z in the form 
([23]) in ((2|) to see what is lacking. We obtain: 

—Dk 2 acos k c x cos q c z + E z (x, z) 
9Ba 2 " 



A + Dk 2 
Ba 3 



16 



a cos k c x cos q c z 



16 



(3 cos 3/c c x cos q c z 



+ cos k c x cos 3q c z + cos k c x cos ?>q c z) 
One can satisfy this equation by putting 



(57) 



III. PROPERTIES OF A SINUSOIDAL DOMAIN 
STRUCTURE 



A + 2Dk 2 c + —Ba 2 = 0, 
16 



(58) 



Now we turn to the discussion of the ferroelectric 
phase. It is quite natural to expect that close to the tran- 
sition the spatial distribution of the polarization may be 
well represented by a single or a few sinusoids with the 



16 A + 2Dk 2 c 
~9 B 



16,4 
~9B' 



(59) 



and introducing the additional "ghost" fields in order to 
compensate the rest of the terms on the right hand side. 
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They are not the real electric fields because they do not 
satisfy the electrostatic equations. In fact, it is inconsis- 
tent to take into account the higher harmonics of the field 
without correcting the polarization self-consistently. But 
all this trouble is not necessary if those "ghost fields" are 
small compared to the real ones. The condition for this 
reads 



Dki > -Ba 



-A< iDki 



1 



:A, 



(60) 



The numerical factor is certainly not reliable and may be 
omitted. 

If the polarization is presented not as a single sinusoid 
but as a sum of a constant and a sinusoidal parts the 
condition of neglecting of the higher order harmonics has 
to be considered anew. In this case, we have close to the 
phase transition: 

P z (x, z) — p + a cos k c x cos q c z. (61) 

For a, we now have from the equation of state, 

= _l 1 A + 2Dk 2 + 3B P 2 

9 B ' y ' 

and among the "ghost" fields there is, for example, the 
component \Bpa 2 cos 2k c x cos 2q c x that can be neglected 
when 



—5 > 7 B P a = B P\ 5 ■ ( 63 ) 



e x k 2 

The maximum of the r.h.s. is at p - 
condition transforms into: 

- A < 2ViDk 2 r . 



B 



-A/6B , so this 



(64) 



which is practically the same as (|60p . Let us emphasize 
that at any temperature the same formula may be applied 
to describe the states with p c < p, since the amplitude a 
cannot be large in this case. 



B. Free energy 

It was convenient and instructive to use the equations 
of state to describe instability of the paraelectric phase. 
However, it is more convenient to use a (non-equilibrium) 
free energy in order to analyze the properties of the sinu- 
soidal domain structures, whose minimization provides 
us with an equilibrium state. The LGD free energy (Q]) 
does not include the energy of the electric field sources, 
and we shall calculate it now. But first, a simple exercise: 
let us perform a similar calculation for a linear dielectric 
with the ideal metallic electrodes at a non-zero voltage U. 



For simplicity, we assume that the dielectric is extremely 
anisotropic and has one component of the polarization 
only, and that only the states with a homogeneous po- 
larization are studied. We have 



F = I dV^P 2 



W Kf 



(65) 



where the integral is taken over the volume of the dielec- 
tric and W e i is the electric energy. The integrand can be 
called the LGD free energy for the linear dielectric. The 
electric energy can be represented as 



W ei 



dV 



El 

8tt' 



(66) 



where the integral is now over "the whole universe" , 
which means in our case over the volume of the dielectric 
and the voltage source. The latter is imagined as a ca- 
pacitor with an infinite capacitance. The electric field in 
the dielectric is always U/l and does not depend on the 
polarization. When we change the polarization (recall 
that we consider the non- equilibrium free energy), the 
charge on the electrodes is changed and, therefore, the 
charge and the energy of the source are changing too. 
Suppose first that the polarization is zero and for this 
state of the dielectric the charge of the source is Q s , so 
its energy is equal to Q 2 /2C, where C is the capacitance 
of the source. Consider now a state with the polarization 
P. The charge on the electrode with the potential U (the 
second electrode has zero potential by assumption) has 
then changed by some value Q. The electric energy of the 



source has changed by 



/2C 



-QU, 



(Q s - Q? - Ql 

where we have taken into account that U = Q s /C and 
that Q 2 /2C — > when C — > oo. In general, the work of 
the voltage source (change of its energy) can be written 
as 

electrodes 

- ^2 e a<Pa = -(eitpi + e 2 ip2) = -e\<fi 

a=l,2 

= -QU, (67) 

where e a are the charges on the electrodes, and we as- 
sumed that the second electrodes in earthed. When the 
polarization is directed outward from the electrode, then 
the charge on the plate is Q = PS, where S is the elec- 
trode (plate) area. Since the electric energy lowers in 
this case, the equilibrium polarization will be directed 
from the electrode U. We have proven a well known fact 
that the polarization tends to be oriented along the field. 
Now, the full change of the energy per unit area due to 
a non-zero polarization is F/S = IAP 2 /2 — PU, where I 
is the thickness of the dielectric plate. Minimizing this 
non-equilibrium free energy with respect to P, we find 
P = 4t- This is, of course, a well-known result for a 
linear dielectric: the equilibrium polarization is directed 
along the electric field and is proportional to the the field 
(cf. the standard P = X-E). Using this, one can find that 
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the sum of the "LGD" energy of a linear dielectric (|65|) 
and the energy of electric field inside it at the equilibrium 



Consider the electric field energy for the sinusoidal 
part. We find from Eqs. and @, that 



A 2 
-P 2 



(68) 



where e = 1 + 47rx, which is another well known formula. 

Now we can evaluate the free energy of a ferroelectric 
capacitor at a given bias voltage. Since we are interested 
in the case of an inhomogeneous ferroelectric polariza- 
tion, we have to take into account the polarization per- 
pendicular to the polar axis, cf. Eq. (|26[). For the LGD 
energy of this polarization together with the energy of 
the perpendicular electric field in the ferroelectric, we 
can use Eq. I|68p . The same is valid for the nonferroelec- 
tric part of polarization of the ferroelectric along axis z, 
as well as for the polarization and the electric field of the 
dead layer. As a result, we obtain a contribution to the 
full free energy: 



E, 



FE 



8tt 



-dV 



FE 



8tt 



dV 



DL 87r 



(69) 



where the first integral is over the volume of the ferro- 
electric (FE) and the second over the volume of the dead 
layer (DL). It is convenient, though not quite consistent, 
to refer to these contributions as to the electric field 
energy 

The total energy consists, therefore, of three parts: the 
LGD contribution, the electric field energy and the en- 
ergy of the voltage source. It can be presented as (here 
and below, we will always imply that the free energy is 
calculated per unit area) 



F 



dV 



FE 



-Pi 



B 



-P 4 



\Da (VxiPz) (VxjPz) 



\ri(d z P z f 
e x E 2 , , e b E 2 



8tt 



8tt 



dV- 



DL 



8tt 



electrodes 

E 



e a fa, 



(70) 



where the first integral is over the volume of the ferroelec- 
tric (FE) and the second over the volume of the dead layer 
(DL). At the moment, we want to study an anisotropic 
film where the polarization changes only along one of the 
transversal axes (e.g. x) and the standard case of q <C k 
(see our discussion above). Then, we can simplify the 
free energy as: 



~ ~ 47r a 

E fz < E fx = — — sinfc c a; sing c z, (72) 

hence, since e b is hardly larger than e±, the energy of the 
electric field in the ferroelectric is 



F 



dV 



FE 



A 2 
— P 2 

2 2 



dV- 



DL 



8tt 



B 4 

— p 4 

4 2 



QU. 



\D{d x P z f 



e b E 2 z , e ± E'; 



8tt 



8tt 



^dV = Aa 2 . 
8ir XJ 2e_i kt 



(73) 



We obtain the field in the dead layer from Eq. (|32|) . Re- 
calling that we are interested in the case kd <§; 1, q <C k 
and q ~ 7r/Z, we find that in the dead layer: 

~ ~ 87T 2 , 4^/2^1/2 

Edx *C &dz — — , , acosk c x — 7-^; acos/c n a;. 



Wd 



e i a 



The energy of the electric field in the dead layer is 



(74) 



Taking into account Eqs. (|43| . (144]) . one sees that the 
ratio of the energies of the electric field in the ferroelectric 
and in the dead layer is n\/3d/d m , i.e. for our case d 3> 
d m the contribution of the deal layer can be neglected. 

Now, it remains to calculate the contribution of the ho- 
mogeneous part of the electric field energy . To this end, 
we use Eq. pT|) for the electric field in the ferroelectric, 
find the electric field in the dead layer from Eq. (fT2"| , and 
obtain: 



8tt fz 



8n~ 



-E: 



e b e e U 2 



dz 



2Trp 2 dl 

le e + e b d &ir (e e l 



e b d) 
. (76) 

The last term here is the unimportant constant. Finally, 
the energy of the external voltage source is —QU and for 
a unit area one finds from Eq. (jllj) : 



and 



Q 



W e ls = — 



I 



I + e b d/e, 



-P- 



I 



-pU = -lpE 



I + e b d/e e 

The total electric energy per unit volume equals 

2np 2 dl 



2c i kl 



le e + e b d 



pE Q . 



(77) 



(78) 



(79) 



Finally, one obtains the LGD contribution to the free 
energy by using Eq. ([BTj) and calculating the integrals for 
(71)9 = f/i- As a result, the free energy per unit area F(p, a) 



takes the form: 



Let us mention that the electric field energies for the 
homogeneous and the sinusoidal parts are additive, as 
follows from the fact that J dx cos kx = . 



A 



£ 2 
—p - 



A 



B 



35 



2„2 



-a p 



9B 
256' 



-pE , 
(80) 
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where 



2Dk 2 c , 



A = A 
£ = And/ (e e l + e b d) 



2Dkj 



(81) 

And/ (ej) - 2Dfc 2 82) 



are the main parameters of the system, Eq = 
e e U/ (ej + ebd) sw U/l the external field for the usual 
case of a thin dead layer e b d -C e e l. The Landau coeffi- 
cients A, B, ... in (fT]) are renormalized by lattice misfit 20 
(see Appendix D), but we do not account for inhomoge- 
neous strain coupled with inhomogeneous polarization in 
Eq. (fTTj) . Its account is irrelevant when one defines the 
point of stability of the paraelectric phase with respect 
to the domain formation. Indeed, the coupling between 
the inhomogeneous polarization and strains is nonlinear, 
while the problem is linear. This coupling begins to play 
a role when the finite amplitude a of the "polarization 
waves" has to be found. Nevertheless, we shall not take 
it into account to avoid unnecessary complications, which 
do not change our main conclusions^ 2 -. 



IV. ELECTRIC PROPERTIES OF SINUSOIDAL 
DOMAIN STRUCTURES 



A. Polarization - External field curves. Different 
critical thicknesses for stability loss and for memory 

By now, we have made an important step of calculating 
the free energy of the film with the dead layers in a single 
harmonic approximation. Already at this level, we can 
study the regions of (meta)stability of the homogeneous 
state with respect to domains and determine the corre- 
sponding conditions. As we have noted a several times 
above, the residual depolarizing field in the ferroelectric 
resulting from an incomplete screening by the electrodes 
tends to split FE into domains, and that defines the prop- 
erties of very thin films. Importantly, not only formation 
but also a response of the domain structure to an exter- 
nal bias voltage appears to be affected by this field by 
a much greater extent than one might have expected^, 
as we will discuss below. For instance, the p — Ej de- 
pendencies have an unusual shape with a negative slope, 
predicted earlier by the continuum theory for a similar 
case at low temperatures^. 

Apparently, the thin ferroelectric films are very well 
suited for an application of the analytical theory of do- 
main structures in ideal crystals. The theory helps to 
understand some experimental data, as we shall demon- 
strate in the next Section. However, at the moment this 
theory is unable to make practical predictions that seem 
most important, like the estimates of the parameters of 
the ferroelectric films that may correspond to an accept- 
able memory performance. This is not surprising, how- 
ever, since there is no consistent theory of the FE switch- 
ing. 

In this regard, it is worth mentioning that the usual 
practice of using the Kolmogorov-Avrami approach to 



describe the switching in FE is very unfortunate. Such 
an approach was designed for considering nucleation and 
growth in the absence of any long range interaction be- 
tween the nuclei, which is far from being the case in 
the ferroelectrics, and without account for the depolariz- 
ing field suppressing creation of isolated nuclei. In other 
words, strong macroscopic Coulomb interaction defining 
the behavior of thin films is completely neglected in that 
approach, and there is no easy way to take it into ac- 
count. On the other hand, considering domain structures 
as a superpositions of several sinusoidal distributions of 
polarization takes consistently into account the Coulomb 
interaction from the very beginning. Such a superpo- 
sition is a readily identifiable periodic domain structure 
for an equilibrium or a metastable state. Considering the 
switching, one most likely should consider a local pertur- 
bation in the periodic structure (nuclei), but it is very 
inconvenient in terms of the delocalized sinusoidal ba- 
sis functions. Nevertheless, discussion of thin films with 
nearly sinusoidal domains gives a possibility to approach 
the problem of switching from a different perspective. 
We shall see that even the first step in this direction, i.e. 
considering purely sinusoidal domain structure with the 
use of Eq. (f80|) leads to the non-trivial conclusions. We 
believe that further description of domain structures in 
this field is promising and highly desirable. 

Here, we present an analysis of the ferroelectric film 
with the sinusoidal domain structure in the single har- 
monic approximation and will carefully indicate where 
this approximation works. By minimizing the corre- 
sponding free energy (|80p with respect to a, one finds 
the equilibrium value of the amplitude: 



«o 



16 , ■ 



3Bp' 



1/2 



4p c nr 



for 



= for \s\ > 1, 



(84) 



where s = p/p c , with p c — y—A/3B the characteris- 
tic polarization. Substituting this solution into Eq.(|80p 
we then arrive at the dimensionless free energy: / = 
3B£ > l~ 1 F depending on the homogeneous polarization 
only through the parameter s: 



f±(s) 



y_ 

3 



v) s2 + tiv 2 ^ - Vv se 



\s\ > 1, 

1*1 < 1, 
(85) 

where e = Eq/<^ is the relative external field, £ = 
£ 3 / 2 /\/3i? is the characteristic electric field, and 



V 



(86) 



is the characteristic temperature (i.e. the relative dis- 
tance of transition temperature from the paraelectric 
phase that depends on the film thickness I). One can 
easily see that the free energy is continuous with first 
derivative with respect to s, while d 2 f_/ds 2 ^ d 2 f + /ds 2 
at s = ±1. 
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The equation of state s 
condition df / ds — reads: 



s (e) obtained from the 



y/y{l~y)s + 



y 3 ' 2 s 3 



y 3 / 2 s 3 



> 1 



M<i, 



(87) 



and the given state is (meta)stable only when d 2 / /ds 2 > 
0. It is very instructive to compare the behavior of the 
FE capacitor with real electrodes (or, equivalently, with 
the dead layers) with that of the system without the dead 
layers, which does show the standard S-shaped polariza- 
tion loop described in the dimensionless units by 



3 /2 

s /y(l + y) S + ^-S 3 = e. 



(89) 



We illustrate the behavior of the equation of state de- 
pending on the parameter y of the system in Fig. We 
have selected different values of y = 1/4, 3/2, and 3, 
corresponding to different temperatures (thus, y = cor- 
responds to the paraelectric-sinusoidal domain structure 
transition for U = occurring at A = 0, at tempera- 
ture Td < T c ), with the results shown in Fig. |4j We 
see that at a relatively small y = 1/4 (for instance, 
for the system not far below Tj) the equation of state 
(|S7p. (155)) has only a trivial solution p = at E = 0, 
i.e. this is the state with no memory. In the field 
e > e c = y/y — 2y 3 / 2 /3 = 5/12 the system is homoge- 
neously polarized, it splits in the lower bias field via the 
second order phase transition into domains with zero net 
polarization in the film at E = 0. At y = 3/2 (Fig. HJd) 
this transition is first order (the change of the order of 
the phase transition occurs at y = 3/8). This is the first 
instance when the state with a spontaneous net polar- 
ization p at Eq = becomes formally possible as a 
solution to the equations of state. However, this state is 
unstable (d 2 f-/ds 2 < at s — ±1), as is evident from a 
negative slope of p — Eq curve for the upper branch of 
/_ in Fig. QJj. 

There is a metastability of a homogeneously polarized 
state (we shall abbreviate it as the MHS) in the re- 
gion 3/2 < y < 3 [d 2 f + /ds 2 > with the free energy 
higher than that of the domain state, /+(±1) > /-(0) at 
Eq = 0] . Formally, both the ferroelectric memory and the 
polarization switching arc possible at these temperatures 
but no conclusion of practical importance can be made 
before calculating the escape time from the metastable 
state. At larger y > 3 (further down from the phase 
transition) the state with the homogeneous polarization 
in the present single-harmonic approximation, Eq. (|6ip . 
has the same or lower free energy than the state with 
p = 0, f+{s m ) < /-(0), where ±s m are the positions of 
the minima of the free energy / + at E = (Fig. O. 
However, this result is approximate. The reason is that 
Eq. (|23p is valid near the phase transition point only. The 
region of validity of this approximation has been esti- 
mated ini£ and above (Sec. IIIA) as roughly —A < Dk 2 , 
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FIG. 4: The equation of state p = p(Eo) for ferroelectric 
film with dead layer in external field Eo for various values 
of the relative temperature y, Eq. I|86|l : (a) y = 1/4, (b) 
y — 3/2, and (c) y — 3. There is a phase transition between 
homogeneous state and the one with sinusoidal domains: (a) 
second order, which becomes first order in cases (b) and (c). 
All those cases are very different from the case of the FE film 
without the dead layer. ( = £ 3y/2 /\/3B is the characteristic 
electric field. 



which means y < 1, at least for the experimental system 
BTO on SRO/STO considered below in Sec. V. 

In a more accurate approximation accounting for 
higher harmonics to describe the inhomogeneous polar- 
ization, the free energy minimum at s = dips lower 
than that of the homogeneous state. It is these higher 
harmonics that convert the sinusoidal domain structure 
into a conventional one with narrow domain walls. For 
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the parts of the curves corresponding to \s\ > 1, these 
higher harmonics are not important (they are when am- 
plitude of the first harmonic becomes substantial) but 
they will change the curves for |s| < 1 substantially. The 
amplitudes of the higher harmonics are to be considered 
as new variational parameters for the free energy, and 
their account will be lowering the estimated free energy. 
Hence, the minimum at p = in Fig.[5]is actually deeper, 
and the homogeneously polarized phase becomes stable 
not at y = 3 but at a larger value (i.e. at a lower temper- 
ature or a larger thickness). Furthermore, it is possible 
that the homogeneous state would always remain less sta- 
ble than the polydomain state in that region. Indeed, in 
the opposite limiting case, i.e. far below the FE transi- 
tion, it has been shown that for any thickness of the dead 
layer the multidomain state has a lower free energy than 
the homogeneously polarized stated. 

Importantly, the energy barrier between the oppositely 
polarized states is strongly reduced by the existence of the 
domain structure at E « 0, Fig. O It is very suggestive 
that one needs to take this into account while considering 
the problem of polarization switching or the problem of 
practical memory, but it is a problem that should be dealt 
separately. We note again that Fig. [4]: highlights that ho- 
mogeneous (single domain) switching is impossible: the 
homogeneously polarized phase loses its stability with re- 
spect to the sinusoidal domain structure at the smaller 
reverse bias (point A) than the reverse bias flipping the 
homogeneous polarization (point B). The two (necessar- 
ily negative) fields can be calculated from Eqs. (|57|) . (|88p. 
One finds the value of e A either Eq.([57|) or Eq. (J55J for 
s = 1: 

2 



->y 



3/2 



y 



1/2 



(90) 



The field cb can be obtained from Eq. (|87| by extrapo- 
lating it to a region s < 1 and calculating s and e corre- 
sponding to the point de/ds = 0. One finds, 

^ I i \3/2 

e B = — (y-i) 



i 



We sec that the difference between the fields vanishes for 
y 3> 1, but it is appreciable for small values of y : thus, 
for y — 2 the difference is about 30%. The applicability 
of Eqs. (jnnj) , dnij) is not limited by the region of the one- 
sinusoidal approximation. 

One may ask a question about what happens in a 
more realistic case where higher powers of P in the Lan- 
dau expansion, P 6 or P 8 , may be important, like in 
BTO on STOi 3 -. In this case the transition seems to be 
second order according to first-principles calculations 23 , 
same as in an earlier parameterization 2 ^, while the re- 
cent parameterization :] 24 ' 25 indicate very weak first order 
phase transition. A full analysis is fairly involved in the 
case of more general LGD functional and is beyond the 
scope of the present work. Nevertheless, the condition 
for domain instability in a simplest approximation takes 
the form 



where F is the LGD functional for the bulk (or a sample 
with ideal metallic electrodes) and the second derivative 
is calculated for p = p e , the homogeneous polarization 
in a given external field Eq (this condition will be mod- 
ified by the strain, but it is reasonable to leave this ef- 
fect out as a first step). Our analysis indicates that the 
metostability of the homogeneous state now starts at a 
smaller y, i.e. closer to the point of the stability loss 
compared with y > 3/2 value for the second order case 
([55)) . However, since a first order transition occurs before 
the stability loss, it is not clear if the temperature inter- 
val between the transition into ferroelectric phase and 
appearance of the MHS increases or decreases. 
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d 2 F/dp 2 = -2Dk 2 c 



(92) 
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FIG. 5: The free energy of the FE film with the dead layer for 
the relative temperature y = 3, when the transition between 
the homogeneously polarized state and sinusoidal domains is 
first order. It is evident that switching proceeds through the 
state with domains with much lower energy barrier for nucle- 
ation (top of the barriers indicated by arrows). An account 
for higher harmonics in Eq. (|61[) will deepen the energy of 
the domain state compared with the one shown here and will 
lower the barrier even further. 



B. Negative slope of P(Ef) curves 

One unusual property accompanying the domain struc- 
ture is the negative slope of the P — P(Ef) curves 
for ferroelectric films with domain structure. Has been 
predicted some time ago for ideal films at low temper- 
atures in Ref. e >>£ and is viewed as a hallmark of a do- 
main structure governed mainly by electrostatics^ 3 .. Ap- 
parently, this was first experimentally confirmed by the 
data from Noh et al*^, who measured the hysteresis loops 
P = P(E ) of ultrathin BTO films on the SRO/STO sub- 
strate. It has been very instructive^ 3 - to replot them as a 
function of a field in the ferroelectric Ef, P — P(Ef). A 
remarkable specific feature of replotted loops is that they 
all have a negative slope, most pronounced at smaller film 
thicknesses. Certainly, the theory does not give hystere- 
sis loops since its application is limited to an equilibrium 
state, but the behavior of the equilibrium curve is clearly 
seen from the data^ 3 -. 
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We shall show below that the negative slope of the 
P = P(Ef) is characteristic of sinusoidal domain struc- 
ture near the phase transition. It is worth noting, how- 
ever, that the negative slope P/Ej < is not an exclu- 
sive property of the domain structure. It already hap- 
pens in the paraelectric phase in the temperature inter- 
val —2Dk 2 c < A < 0, where the paraphase is still stable 
with regards to small inhomogeneous perturbations. To 
see this for the sinusoidal domain structure, we express 
the external field through the field in the ferroelectric Ef 
with the use of (jlip the additional dimcnsionlcss param- 
eter ijj = 1 + £ (e e l + e^rf) /And = nd m /\/3d. We obtain: 



Andp & 

Hjf — £yQ — ; — &Q 



ej + e b d 



1-ip 



(93) 



To find the response of the sinusoidal domains to the field 
in the ferroelectric, dp/dEf, it is most illustrative to use 
the equation of state analogous to (|88|) written in natural 
units, like 



.4 



£ \p-3Bp 3 



\p\ < Pc 



(94) 



Now, it is convenient to divide this equation by Bp\ = 
—Ap c and express E through Ef from Eq. (f9"3")l . We can 
then rewrite (l88l) as 



1 _ ip \ P_ (p_ 

3 V(l-ip)) Pc \Pc 

This yields for the slope at Ef = : 
dp 1 3y(l - ip) l-ip 



E 



f 



3Bp 3 c 



(95) 



dEf Avii-ty-W C y(i -$)- 3^ 

3eJ 1 



4ndy(l-ip) - 3ip' 



(96) 



where we have used y — —A/£, and £/(l — ip) = 
4ird/(e e l + ebd). Finally, this gives for the "dielectric con- 
stant" of the ferroelectric film : 



e f = 1 + AndP/dEf 



3e e l 



1 



d -V) -3^' 



(97) 



At the point of stability loss with respect to sinusoidal 
structure, the dielectric constant is already negative, as 
we mentioned above: 
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This result is, of course, identical to the standard value 
of the "dielectric constant" of the FE film itself 
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(99) 



at the border of stability loss into sinusoidal domain 
structure. Note that the negative sign of e/ does not 



mean an instability, because this coefficient does not 
characterize the response of the system as a whole to 
an external field, which is Eq, not Ef. The capacitance 
of the device is defined as a response to Eq that is always 
positive^. 

One can see that the slope becomes positive at y > 
3ip/(l — ip), yet this result is beyond domain of appli- 
cability of the single- harmonic approximation. In fact, 
we know that the slope is negative for the film with do- 
mains at low temperatures, as has been predicted some 
time ago^: it is a hallmark of domain structure governed 
mainly by electrostatics. Indeed, a net polarization of 
the domain structure is P = in zero external field 
Eq. At small Eq there will be positive net polarization 
in external field because of growth of domains oriented 
along the external field, and the resulting negative field 
in the FE. Thus, the "dielectric function" of the film is 
negative, e/ = 1 + AndP /dEf\E f =o < Ot Comparing 
theory with the data, the expression for the "dielectric 
constant", given by Eq. (31) of RcfJ., can be simplified 
to 



-e e l/d. 



(100) 



Substituting the numbers for the BTO on SRO/STO film 
with thickness / = 5nm, we find the theoretical value £/ = 
—525 for equilibrium conditions while the experimental 
one found from Fig. [6]lb (raw data for 2kHz, Ref.— ) is 
e/ = —680, i.e. both estimated and measured values are 
pretty close. The slope is negative in all samples with 
thicknesses up to 30nm studied in Ref^, see Fig. [71 



V. EXPERIMENTAL EXAMPLE: BATIO3 
FILMS WITH SRRUO3 ELECTRODES ON 
SRTIO3 SUBSTRATE 

Our above assumption of a second order phase tran- 
sition may be a reasonable approximation for BaTi03 
films on SrRu0 3 /SrTi0 3 (BTO/STO) substrate that 
have been investigated down to thickness / = 5n m 14 ' 15 ' 16 . 
We use it to illustrate our theoretical results while dis- 
cussing, in particular, the effects of the higher order 
terms in the LGD free energy. Because of —2.2% lat- 
tice mismatch, the BTO film effectively becomes a uniax- 
ial FE and can be treated within the Landau-Ginzburg- 
Devonshire (LGD) free energy very accurately (results 
are precisely the same as in ab-initio calculations where 
the latter are available^.) 

Our aim is to describe what we think to be a consis- 
tent method of interpreting the experimental data, i.e. 
the ferroelectric hysteresis loops for the BaTi03 films 
on SrRuOa/SrTiOs with different thicknesses, from 5 
to 30nm. The electrode parameters were found to be^: 
d/2 = A = 0.8A, e e — 8.45. The parameters of the mate- 
rial, i.e. the LGD coefficients for the homogeneous states 
are taken from^ with a renormalization due to the lattice 
misfit strain and the clamping according to2£, see the Ap- 
pendix A. This set is different from the one used in Refi^ 
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FIG. 6: (a) The external field E ob (where E f = 0) and 
Ear where the relaxation of polarization starts in 5 nm thick 
filrn^. Inset: the spontaneous polarization P s and the ex- 
trapolated P(r^. (b) The measured P(Eo) and the "actual" 
P(Ef) hysteresis loops. 

and there are qualitative differences in some experimen- 
tally verifiable conclusions for the two sets. In particular, 
the coefficients of^ imply a weak first order paraelectric- 
ferroelectric phase transition, while those of^S imply a 
second order one, which is also suggested by the results 
of the first-principles modeling^. Since the definitive ex- 
perimental studies were not carried out, it is not clear 
which set is better, and we have chosen to use the set of 
Refi^, when discussing homogeneously polarized state, 
but the assuming second order phase transition while 
treating the domain structures. 

What remains undefined are the additional boundary 
conditions (ABC) parameters, but we show below that 
they can be neglected in the case of BTO on SRO/STO. 
This can be seen from the analysis of the experimental 
data for single domain states. The role of the ABC for 
homogeneously polarized states has been studied in many 
papers in the present geometry when the polar axis is per- 
pendicular to the film plane, starting with Kretschmer 
and Binder—. Specific features of a ferroelectric sur- 
face/interface have been taken fully into account only 
recently-^. For a symmetric system (two identical sur- 
faces) and single domain state their effect is reduced to 
a renormalization of the coefficient A, similar to that of 
Ref J£ A -> A\ = A + (2a + (3) /I, where a and (3 are the 
parameters pertaining to a surface/interface. One can 
find the polarization directly from the experiment which 
would correspond to a zero external field, Eq = 0, if there 
were no domains. In Ref it has been found by extrap- 
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FIG. 7: The measured P(E ) and the "actual" P(E f ) hys- 
teresis loops: (a) film thickness I = 9 nm, (b) I =30 nm. 



olation of the high external field part the p (Eq) curves 
to Eq = 0. Generalizing Eq. (TP5|) to account for the ABC 
and the higher order terms in the LGD free energy, we 
obtain: 

(Ax + P + Bp" + Cp 5 + Fp 7 = Eq, (101) 

i.e. for Eq = 

Ax + ^ + Bp% + Cp% + Fpl = 0, (102) 

and one can find Ax (I) from the experimental data for 
Po (I). The equation for spontaneous polarization P s (po- 
larization in the film with zero internal field, Ef = 0): 

Ax + Bp 2 s + Cp 4 s + Fpl = 0, (103) 

and find it as a function of /. Since wc know from the data 
the polarization at zero external field, po = p(Eq = 0)i^, 
we can easily obtain p s from (|102p , (|103|) . The functions 
p e o (I) and p s (I) are shown in Figd^inset). We see that 
the dependence p s (I) is weak, i.e. the ABC can be ig- 
nored, since the variation of the spontaneous polarization 
P s with thickness due to the (2a + f3)/l term is small, and 
we shall assume Ax — A. 

Let us now discuss what the theory tells about the state 
of films of different thickness at Eq = 0. Let us consider 
first the case of B > 0. This case has been discussed 
in Sec. IV. We shall start with this case and then try 
to generalize it. From Eq. (|10ip for Eq = 0, one sees 
that the solution p = (the paraelectric phase) becomes 
the only possibility at I < And/ (— Ae e ) . Using the above 
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value for the coefficient A, we find that it happens at 
I = i/j = 3.5nm. 

However, we also need to study the stability of the 
paraclcctric and homogeneously polarized state with re- 
spect to the domain formation. First, we calculate the 
value of critical dead layer d m for our system. The value 
of the transversal dielectric constant, ej_, can be calcu- 
lated with the use of the LGD coefficients, and we find 
that ej_ = 218 at room temperature (Appendix D). Then, 
using Eq. we find that d m — 0.34A, i.e. in our 

case d = 5d m . Using values of A, e e from Refj^ 4 -, calculat- 
ing e± using the coefficients of Ref. 24 (Appendix D), and 
the value of ^/D/An = 0.2A from Ref. 26 (Appendix E), 
we obtain the following critical thickness for domains at 
room temperature (ej_ = 218): 



RT 



(f 



1/2 



/\A\ ~ 1.6 nm. 



(104) 



Since d > d m , then ^ T < Z^ T , and the phase transi- 
tion is into a multidomain state. The spatial distribution 
of a spontaneous polarization is near sinusoidal at I > 
lj T . Higher harmonics develop with increasing thickness, 
and the polarization distribution tends to a conventional 
structure with narrow domain walls. The half-period 
of the sinusoidal domain structure can be estimated as 
a r = 1.7 nm at the transition and as a RT = 2.2nm for 



I = 5nm from a ss Tr/k c and (|45[) . 

It is instructive to consider the phase transition with 
regards to the thickness at zero temperature, where we 
get e± = 408, ^ 0K) = 0.8nm, 1 { ° K) = 2.5nm. The last 
result (a homogeneous critical thickness of 2.5nm) is re- 
markable, since it practically coincides with the ab-initio 
calculation for the critical thickness of 2.4nm in Ref. 23 . 
The ground state of the film is, however, not homoge- 
neous but multidomain, and the domain ferroelectricity 

appears in films thicker than l^ K ^ = 0.8nm, which is the 
true "critical size" for ferroelectricity in FE films in the 
present study at zero temperature. 

To use the results of Sec. IV, it is convenient to ex- 
press the variable y as a function of the film thick- 
ness. Since 2Dk 2 oc one can put 2Dk 2 (I) = 

l c 2Dk 2 c {l c ) jl = -Al c /l. Then A (I) = A - Al c /l = 
A i__L^ xhe parameter £ (/) = 4^/ + e&c Q _ 2Dk 2 (I) 

can be obtained as follows. The phase transition into 
the homogeneous state would occur at I = at 
A + 4nd/ (ej h + e b d) = (Sec. IIB). Therefore, A = 
—4ird/e e lh, hence 47rd/ (e e l + e^d) ~ iird/ej — —Alh/l 
and £ = -Alh/l + Al c /l = -A(l h - l c ) jl. Finally, we 
get y = {I — lc) I (Ift ~ lc)- We see that the beginning of 
the memory for the case of a second order transition from 
the paraelectric phase, which is at y = 3/2, corresponds 
to I = 4.5nm. Recall that for this case the homoge- 
neous phase becomes stable at y > 3 (or for I > 7.3nm) 
within our approximation, but this conclusion is unreli- 
able because, as we discussed above, the homogeneously 
polarized state will almost certainly remain metastable. 
Returning now to the experimental data for BTO on 



SRO/STO, we should mention that the thinnest film 
(5nm) is thicker than the academic "limit for memory" 
even if B > and even more so for B < 0. This does 
not mean, however, that the sample has a real stable 
memory because apart from the memory loss because of 
absolute instability of the homogeneously polarized state 
there is also a memory loss because of domain nucle- 
ation. There is still no reliable theory of this process, but 
we can gain some insight into the process by analyzing 
the experimental data. No memory has been found with 
a lifetime longer than 1 second in BTO on SRO/STO 
capacitors with thicknesses less than 30nm 14 . Kim et al 
have found that during the observation time t re i ax = 10 s 
the domains begin to form only if the applied external 
field along the polarization is less than Eq t - — 910kV/cm 
for 5nm sample (raw data). We can immediately esti- 
mate from the results above that at Eq = E$ r the do- 
mains begin to form only when there is a field in the film 
Ef = — (490 ± 70) kV/cm opposite to the polarization. 
This is natural: it was known for a long time that in or- 
der to nucleate domains one has to apply an "activation 
field" (E a )21. Studying relatively thick BaTi0 3 samples 
with I = (2.5 - 35) x 10 _3 cm, Mera 2 ! found that this field 
depends both on the waiting time r and on the sample 
thickness I and have obtained an empirical formula for 
it: 



E a = (alnr) jl. 



(105) 



In his thick samples, Merz identified the activation fields 
to be up to 20 kV/cm, which looks quite high if one were 
to apply the Merz's scaling law (|105[) to Noh's samples 14 . 
Indeed, the thicknesses of the films in Refj^ 4 - are four or- 
ders of magnitude smaller, while the fields are less than 
the two orders of magnitude larger. It is evident that 
the Merz's formula cannot be applied literally to the 
films of all thicknesses. (To be fair, Merz studied not 
the strained BaTiOs, but this hardly accounts for such 
a difference). However, within the interval of thicknesses 
studied in Refi 1 ^ the Merz's formula describes, at least 
qualitatively, the dependence of the activation field on 
the film thickness (see FiglH]). It is worth mentioning 
that at small thicknesses the domain structure is nearly 
sinusoidal, and one can expect weaker pinning compared 
to thicker films. It is hardly surprising that the empirical 
Merz's formula obtained for conventional domain struc- 
ture does not apply to a sinusoidal one. 

Using Eq. (jlOip . one can find that the electric field in 
short-circuited sample is about 1200 kV/cm, exceeding 
the magnitude of the estimated activation field. This 
means that in a short-circuited sample single domain 
state relaxes faster (perhaps much faster) than in 10 3 s. If 
the value of the activation field is defined by the thickness 
only and not by properties of electrodes or an electrode- 
film interface, one can speculate about the properties of 
electrodes that can facilitate a smaller field in a short- 
circuited sample and a longer, at least 10 3 s, retention 
of a SD state. We have found that for I = 5nm such 
an electrode should have X/e e = d/2e e < 0.043A. In- 
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deed, with such a hypothetical electrode the depolarizing 
field in an unbiased 5nm sample would be less than 500 
kV/cm discussed above. Since in the BTO samples^ the 
value A/e e is about O.lA, it does not seem totally impos- 
sible to find such an electrode. This will correspond to 
d = 0.73A« 2d m , the phase transition would be still into 
a multidomain state. The homogeneously polarized state 
that is metastable at lh = 1.4nm would be retained for 
at least 10 3 s in the 5nm BTO film with such an electrode 
according to the discussion above. One can estimate the 
critical thicknesses for the domain states to see that here 
too Id < lh (since d > d m ), although the estimates should 
use more accurate formula for k c , compared to (|43|) . be- 
cause d is getting close to d m . While looking for possibil- 
ities for ferroelectric memories, it seems more promising 
to optimize both the properties of the electrode and the 
FE material. Recall that d m oc D 1 ^ 2 and the constant D 
proved to be particularly low for the (100) orientation of 
BaTiO^. From a theory viewpoint, we see that one can 
play with electrodes, materials, orientations of the film 
growth, and the misfit strains to find a combination that 
may be suitable for the nanosize ferroelectric memory 
elements. 



VI. CONCLUSIONS 

We have described the application of the Landau the- 
ory to a problem of (meta)stability of a homogeneous 
state in thin ferroelectric films. The theory allows to bet- 
ter understand the data and gain more insight into the 
practically important problem of memory (meta)stability 
in thin films. One advantage is that the theory can start 
with relatively simple ground state. We have proven, 
in particular, that switching of a polarization neces- 
sarily goes via the intermediate domain state, the em- 
pirical fact well known from experiment. The emerg- 
ing domain structures in thin films are quite unusual, 
the "sinusoidal" ones. We show the distribution of the 
polarization, which is reminiscent of the one for do- 
mains with abrupt distribution of polarization (the Kit- 
tel structure^.) To slow down the polarization relaxation 
into a domain state, one can select different electrodes or 
ferroelectrics, and the theory guides one as to what pa- 
rameters are favorable. For instance, ferroelectrics with 
large energy of domain walls (large gradient, or inho- 
mogeneous terms in the free energy) should have better 
memory retention. 

There are many open questions for both theory and 
experiment. For instance, there is no acceptable theory of 
the polarization switching as of yet and, therefore, theory 
cannot predict what parameters are needed to achieve a 
needed memory retention. How the sinusoidal domains 
evolve with temperature and/or thickness of the film into 
the standard domains with sharp domain walls, remains 
to be investigated. Obviously, the continuous theory will 
play very important role in this future work. 
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APPENDIX A: HISTORICAL NOTE 

The mathematical equivalence of electrostatics and 
magnetostatics led, at the beginning of development of 
the theory of ferroelectric domain structure, to apply- 
ing the relevant results of the theory of magnetic domain 
structures to the case of ferroelectrics. The ferromag- 
netic domains were postulated by Weiss in 1907. How- 
ever, until 1935 their origin had not been actually under- 
stood. Several authors tried to understand their origin 
within statistical mechanics of an infinite medium. This 
was even the case with such an outstanding scientist as 
F.Bloch^ curiously, in the same paper where he consid- 
ered the famous "Bloch domain wall". In 1935 Landau 
and Lifshitz 1 pointed out futility of these efforts and indi- 
cated the demagnetizing effect of surfaces as the origin of 
ferromagnetic domains, i.e. the demagnetizing field ex- 
isting in finite magnetic bodies. Landau and Lifshitz pro- 
posed a fairly complete theory of ferromagnetic domain 
structure and obtained some formulas which since then 
were often attributed to other authors. For instance, the 
proportionality of the width of the domains to the square 
root of a film thickness I is frequently attributed to Kit- 
tel (see alsc4). Kittel highly appreciated the Landau and 
Lifshitz work and has further developed their theory in 
1940's summarizing the Landau and Lifshitz's and his 
own results in his often cited review^. For later devel- 
opment of the ferromagnetic domain theory see, e.g., a 
textbook 3 ^. 

In the theoretical part of their pioneering paper on the 
domain structure of Rochelle salt and KH2PO4, Mitsui 
and Furuichi 3 pointed out the relevance of the earlier re- 
sults for ferromagnets with strong uniaxial anisotropy to 
their case. Also in 1950's the group of Kanzig 3 - clearly re- 
alized that the depolarizing field in ferroelectrics should 
be much more important than in ferromagnets. They 
argued that in a slab of an uniaxial ferroelectric, e.g. 
KH2PO4, with the polar axis perpendicular to the slab 
plane no phase transition into mono-domain state is pos- 
sible and the only possibility of phase transition into fer- 
roelectric phase is with formation of a domain structure. 
This phase transition should occur at a lower tempera- 
ture and the lowering of temperature is inversely propor- 
tional to the slab thickness. They tried to estimate this 
lowering of the phase transition temperature to find what 
is now called "the critical thickness for the ferroelectric- 
ity" in the case of KH2PO4. They also clearly under- 
stood that their estimates are irrelevant to the case of a 
three-axial ferroelectric such as BaTi03 where the closure 
domains (nowadays called "vortices") should form. The 
possibility of screening of the depolarizing field not only 
by free charges from the electrodes but also by the charge 
carriers of the material itself was also clearly realized at 
the very beginning of a systematic study of ferroelectric- 
ity One can see the relevant discussion, in e.g. papers 
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by the group of Kanzig 31 . 

Ivanchik 32 considered the possibility of a mono-domain 
state in a slab of nonelectroded uniaxial ferroelectric tak- 
ing into account the screening of the depolarizing field by 
the carriers of the material, which was treated as a semi- 
conductor with a large band gap. Later he and others^ 3 - 
took also into account the screening by metallic and semi- 
conductor electrodes. Also in 1960s, the effect of finite 
screening length of metallic electrodes on dielectric mea- 
surements has been discussed by several author o 34 ' 35 ' 36 . 
Batra and Silverman 3 ^ criticized some details of an ap- 
proach adopted by Ivanchik's group and considered an 
effect of incomplete screening by the electrodes on the 
phase transition into the ferroelectric state assuming that 
the phase transition proceeds into a mono-domain state. 
This topic was later developed in a series of papers by 
Batra et aZ.— . Similar to the Ivanchik's group, the main 
assumption of two series of paper, stability of the mono- 
domain state was never checked although earlier Bjork- 
stam and Oettel 39 pointed out that incomplete screening 
due to a near electrode dead layer may lead to the same 
domain structure as in non-electroded samples. It is clear 
now that the assumption about a monodomain state is 
almost never satisfied for realistic parameters of the sys- 
tems, at least close to the phase transition. Chensky 40 
substantially improved the Kanzig's group determination 
of lowering of a phase transition temperature due to im- 
possibility of transition into a monodomain state in a 
nonelectroded slab of a ferroelectric with the polar axis 
perpendicular to the slab. Selyuk considered screening of 
the depolarizing field by both space and surface charges 
and found that below the phase transition a monodomain 
state may become energetically more favorable than the 
multi-domain one 41 -. 

The next important step was the Chensky and 
Tarasenko paper mentioned in the Introduction. Re- 
cently, the Chensky- Tarasenko approach to study dif- 
ferent cases of the domain structure formation at phase 
transitions has been used by Bratkovsky and Levanyuk 8 
and by Stefanovich et al. i2 for the case of phase tran- 
sition in ferroelectric periodic structures. Some other 
papers are cited in the main text. The depolarization 
field screening due to charge carriers of the material and 
those of the electrodes has been recently reconsidered 
by Watanabe 43 . He obtained results substantially dif- 
ferent from those of the Ivanchik's group at the expense 
of making several important mistakes, he also avoided 
any comparison of his results with those of the previous 
authors. 

This historical note does not pretend to represent an 
exhaustive review of development of theory of the do- 
main structures in ferroelectrics. We have not mentioned 
many important works aiming at citing only the most 
important ones falling into the scope of the present pa- 
per. In particular, we did not mention many papers de- 
voted to domain structures far from the phase transition 
(one could call this case the "Kittel limit" versus the 
" Chensky- Tarasenko limit" ) as well as papers related to 



fervoelastic domains. Although conceptually the latter 
topic is close to the theory of 180-degree domains forma- 
tion considered in this paper, it is still well beyond its 
scope. 

APPENDIX B: ADDITIONAL BOUNDARY 
CONDITIONS 

The equation of state ([2|) relates electric field at a given 
point not only to the polarization at the same point but 
also to its spatial derivatives. In the electrodynamics of 
continuous media this is referred to as an account for 
spatial dispersion of dielectric constant 4 -,—. Differential 
equations of the electrodynamics with account for the 
spatial dispersion are of higher order than in usual ("lo- 
cal") electrodynamics, i.e. to specify among solutions 
of these equation those that correspond to physical sit- 
uation one needs more boundary conditions than in the 
local electrodynamics. These are the "additional bound- 
ary conditions" (ABC). Their origin for different physical 
systems is discussed in detail in^ 4 -. 

Another way to ABC consists in applying the Landau 
theory of phase transitions not to an infinite medium 
but to finite systems. For the first time such an ap- 
proach has been proposed in the paper by Ginzburg and 
Landau 4 ^, where the Landau theory has been used, in 
particular, to analyze properties of thin superconduct- 
ing films. The authors argued that at the boundary 
superconductor-vacuum or superconductor-dielectric the 
derivative of the order parameter along the normal to the 
surface should be zero in the absence of a magnetic field. 
With such boundary conditions there is no dependence 
of the phase transition temperature on the film thick- 
ness. As far as we know, the first example of the phase 
transition temperature dependent on the film thickness 
was given by Ginzburg and Pitaevsky 4 ^, who considered 
superfluid phase transition in thin films of liquid helium 
sandwiched between two solids. They argued that the 
order parameter should be zero at liquid helium-solid in- 
terface. The influence of this boundary condition can be 
easily understood considering the loss of stability of the 
normal phase in the same way as in the present paper. 
Replacing in Eq. (fT]) P z by -q one obtains the free energy 
of the Landau theory of phase transitions 4 ^. Since in the 
case of superfluid transition there is no physical field con- 
jugated to the order parameter instead of Eq. (fT7j) one 
has 

Ar, - gdlr) = 0, (Bl) 

with the boundary conditions r\ (z = ±Z/2) = 0, where I 
is the film thickness. Since at A > the solutions of this 
equation are real exponentials, no non-trivial solution 
satisfies the boundary conditions and the normal phase 

is stable. At A < 0, it has a solution cos (^/ —A/ gz^j 

which satisfies the boundary conditions, if y^—A/gl = it. 
This value of A corresponds to the loss of stability of 
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the normal phase, i.e. the phase transition takes place 
at A(T = T cf ) = -it 2 'g- 1 ^ 2 . Within the Landau the- 
ory A — A' (T — T c b) , where T c b is the bulk transition 
temperature. Therefore, T c f — T c (, cx — l~ 2 . 

In 1960s, there was an intensive discussion of the 
"proximity effects", i.e. the properties of a system con- 
sisting of a normal metal deposited on top of a super- 
conducting metal. For the superconductor film the main 
effect was the lowering of the phase transition temper- 
ature. At the end, it has been realized that for suffi- 
ciently thick layers of the two materials and in the ab- 
sence of magnetic field this effect can be described within 
the Ginzburg-Landau theory, i.e. the Landau theory for 
superconductors completed by the boundary conditions: 

Ag±. = 0, (B2) 

for z — ±1/2, where I is the thickness of the superconduc- 
tor film sandwiched between layers of a normal metal and 
A is the parameter of the interface, the so-called extrap- 
olation length. Boundary conditions of this form were 
proposed by de Gennes ? while considering another sys- 
tem, the Josephson junction, while their general char- 
acter was understood by Zaitsev^. These de Gennes- 
Zaitsev boundary conditions^ proved to be applicable 
well beyond the theory of superconductivity and were 
derived later on for many different systems. 

For A > 0, the only case considered by de Gennes and 
Zaitsev, the phase transition in the film occurs again at a 
lower temperature than in the bulk. Indeed, the solution 
of Eq. (jBlj) is again cos {^J —A/ gz^j but to satisfy the 

boundary conditions the value of A should satisfy the 
equation 

tan (ki//2) = — !— , 
Aki 

where k\ = y/—A/g. For / ^> A, the solution of this 

equation is A ~ — i^g^ 1 (I + A) 2 and, therefore, the 
form of distribution of the order parameter at the phase 
transition point is cos [irz/ (I + A)] , and this functional 
form explains the term " extrapolation length" . 

Later, these conditions appeared in study of ferro- 
magnets derived from purely phenomenological point of 
view by Kaganov and Omelyanchuk 51 . These authors 
have taken into account dependence of the surface en- 
ergy of the crystal on the order parameter via terms 
arj 2 (±2/2) /2 in the free energy per unit surface consid- 
ering ferromagnetic phase transition in a slab with the 
anisotropy axis in the slab plane. In this geometry, no 
magnetic field arises due to inhomogeneous magnetiza- 
tion. As a result, they have arrived at the condition 

ff ^±a?7 = 0, (B3) 

which is the same condition as Eq. (|B2[) . The most im- 
portant novelty of Refi^i is that a possibility of the two 



signs of A = g/a ("positive and negative surface energy") 
has been discussed. In the case of negative A, a non- 
trivial solution is possible at A > 0. It is cosh 'A/ gz^j , 
which satisfies the boundary conditions if 

tanh(K2'/2) = t — i 

where n 2 = y/A/g. For I ^> |A| , this equation has the 
solution K2 = |Ap , i.e. at A = g\ 2 . We need to men- 
tion that now the phase transition occurs at A > and 
the temperature of the transition is independent of /, i.e. 
the same transition would occur at the surface of an in- 
finite half-space. The form of distribution of the order 
parameter at the transition is characterized by exponen- 
tial fall-off of the order parameter while going from the 
surface into the bulk, i.e. the phase transition is now 
not a bulk but a surface transition. The width of the 
region affected by the transition increases approaching 
the temperature of the transition in the bulk, so that all 
the results for the film are to be affected by the tran- 
sition. For an infinite half-space, this occurs at A = 
only. Of course, the Kaganov-Omelyanchuk treatment is 
valid also for ferroelectrics with the polar axis lying in 
the slab plane as well as for any order parameter which 
is not coupled with long-range fields. 

Kretschmer and Binder 18 , while using the same bound- 
ary conditions, have taken into account that for ferro- 
magnets and ferroelectrics with the magnetization (po- 
larization) perpendicular to the surface the surface ef- 
fect will be affected by magnetic (electric) field arising 
due to the inhomogencities in the magnetization (polar- 
ization). For ferroelectrics, their case is that of a slab 
with short-circuited ideal metallic electrodes, i.e. of a 
transition into single-domain state. They found essential 
differences with the case of absence of the depolarizing 
field, i.e. they found essential differences between the 
cases of polarization parallel and perpendicular to the 
surface. Tagantsev and Guerra^ have demonstrated re- 
cently that in this problem one needs to take into account 
the noncritical part of the polarization. We shall follow 
them in this regard. 

Along with a modified Eq. (|Bip . which has now the 
form 

AP - gd 2 P = E, (B4) 
one has to consider an electrostatics equation 
d z {e b E + 4ttP) = 0, 

or 

47T 47T f l/2 

E + —P = C=— Pdz, (B5) 
£fc e&t J -1/2 

where the constant C is found from the condition of the 
short circuit: J_^ 2 Edz — 0. Substituting Eq. (|B5[) into 
(|B4[) one obtains 

{A + 47r/e h ) P - gd 2 P = C, 
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and has to look for a non-trivial solution of this equation 
satisfying the boundary conditions (|B2|) . The form of the 
solution does not depend now on the sign of A, it is 



P 



C 



A + 47r/e b 



C\ cosIikz, 



where k = y/(A + 47r/e{,) Jg and C\ is another constant 
related to C via Eq. (IB5|) . As a result, one finds 



P = 



C 



A + 4ir/e b 



AkI 



2sinh(^/2) 



cosh/ 



The boundary conditions (|B2|) are satisfied if 



.4 



kI (Xk + coth(«//2)) 
2a 



kI (Xk + 1) 



kI (g\K + a) 



(B6) 



Since K _1 is the same order of magnitude as the inter- 
atomic distances, the approximation in this equation is 
very good. The polarization distribution at the transition 



P <x 1 



(gn + a) sinh (kI/2) 



cosIikz. 



(B7) 



One sees that now the phase transition occurs in the 
whole volume, both for a "hard", a > 0, and for a "soft" , 
a < 0, with the polarization that is practically constant 
over the volume being somewhat smaller for "hard" sur- 
faces and somewhat larger for the soft ones. The shift 
of the phase transition temperature in the both cases is 
proportional to 

Since K _1 is very small, the formulas for the polariza- 
tion distribution near ferroelectric surfaces are beyond 
the continuous medium theory used by Kretschmer and 
Binder and should be considered as qualitative. In this 
sense account for the noncritical polarization is not very 
essential here. Another shortcoming of their theory men- 
tioned above is overlooking the "polar" character of a 
surface which is not important for ferromagnets but is 
essential for ferroelectrics. This importance has been 
emphasized by Levanyuk and Minyukov 53 . While dis- 
cussing surface structural phase transitions, these au- 
thors have pointed out that when rj can be identified with 
polarization perpendicular to the surface, the de Gennes- 
Zaitsev boundary conditions should be complemented by 
an additional term which makes them inhomogeneous. 
Bratkovsky and Levanyuk^ have explicitly taken into 
account the complete boundary conditions considering 
phase transition in a slab of uniaxial ferroelectric with 
ideal metallic electrodes. Because of the inhomogene- 
ity of the boundary conditions, the polarization is never 
zero over all the volume, even in the paraelectric phase, 
though the region where it is not zero is the same as the 
near-surface region in the Kretschmer and Binder theory, 
i.e. its thickness is comparable to interatomic distances. 
Still, to study stability of the paraelectric phase one has 



to linearize the problem not with respect to P = 0, but 
with respect to P e (z) , an inhomogeneous distribution 
existing in the paraelectric phase. As a result, instead of 
Eq. (|B4[) for P one has an equation for 5P = P — P e (z) 
with the coefficient of the first term A + BP^ (z) instead 
of A in Eq. (|B4|) . This means an additional "hardening" 
of the surface layer or an effective positive renormaliza- 
tion of the coefficient a calculated in Ref.-S. This means, 
e.g., that a surface which is "soft", a < 0, for polariza- 
tion parallel to the surface, may be "hard" with respect 
to polarization perpendicular to the same surface. The 
polarity of surfaces leads also to smearing of the phase 
transition if the two surfaces are different^. 

Here, we have given only a very brief overview of the 
phenomenological theory of the surface effect on ferro- 
electric phase transitions leaving aside both microscopic 
theories and discussion of other effects of the surfaces, 
e.g. their role in the polarization switching (see^). In 
1990s, when thin ferroelectric films became a popular 
topic there appeared numerous papers reconsidering the 
phenomenological theory and trying to find its parame- 
ters from the experimental data or microscopic theories. 
We shall not review these papers because, to our under- 
standing, they did not advanced much the development 
of the phenomenological theory in comparison with the 
papers cited above. This certainly true of the formula- 
tion of the boundary conditions and finding the phase 
transition temperature that we were concerned about in 
this Appendix. 

We do not take into account the ABC in the present 
review. From the point of view of electrodynamics of 
continuous media, this is allowed since the "spatial dis- 
persion of the dielectric constant" is taken into account 
in the direction parallel to the surface (x) and, in ef- 
fect, not in the direction perpendicular to the surface 
(z): the coefficient g does not enter our final formulas. 
Physically, this is justified by the rigidity of polarization 
with respect to its variation along z— axis, which we have 
discussed while reviewing the results of Kretschmer and 
Binder. In addition, the experimental data for the sys- 
tem (BTO on STO/SRO) we have used to illustrate the 
theory indicates that the ABC can be neglected^. 



APPENDIX C: METHODOLOGICAL NOTE 

There is a lot of confusion in the literature with regards 
to the formula for the energy of a ferroelectric with ac- 
count for both the voltage source and the electric field 
due to the polarization of the material. That is why we 
preferred to use the equations of state as long as possi- 
ble and, when being forced to consider the free energy, 
we demonstrate its correctness by comparing the results 
obtained from the equations of state and from minimiza- 
tion of the free energy. By the way, the trouble with the 
Watanabe's papers^ 3 - mentioned in the Historical note is 
due to his use of the free energy while the Ivanchik's 
group used the equations of state to treat the same prob- 
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lem. 

To begin with, we note that some authors write the 
LGD functional as a function of P while others did it as 
function of the dielectric displacement D. The first group 
includes Ginzburg^, Devonshire^, Jona and Shirane^, 
as well as many others, e.g. Batra et al. and the present 
authors. The second group includes importantly Lines 
and Glass^ and, e.g. Ivanchik et al^^. There is also 
the third group led by Tagantse v 59 ; 52 who pointed out 
that in some situations it is important to take into ac- 
count that the order parameter does not refer to the to- 
tal polarization but only to a part of it and there is also 
another "background" or " noncritical" contribution. In 
some problems it should be taken into account when the 
effects of depolarizing field are considered. In this review, 
unlike in our previous works, we have taken it into ac- 
count but it proved out that it disappears from the final 
results for the problems that we considered here, which 
is not necessarily the case for some other problems. 

We think that the origin of existence of these variations 
in the functional representations is the Landau and Lif- 
shitz book " Electrodynamics of Continuous Media" pub- 
lished in the USSR in 1957 and translated into English 
in 1960^. This excellent book had a mishap of using 
D instead of P exactly in the section devoted to ferro- 
electrics. It has been eliminated in the second edition of 
the booki, at least partially, but, because, perhaps, of a 
great authority of Landau, the confusion originated. 

The reason why Landau and Lifshitz put what we call 
now the LGD thermodynamic potential as a function of 
D is that the equilibrium thermodynamic properties of 
a system in a presence of an external voltage source are 
described by thermodynamic potential depending on D. 
This reasoning is flawed, however, because the thermody- 
namic potential (or free energy) of the Landau theory of 
phase transitions^ is a non- equilibrium thermodynamic 
potential (free energy) . It depends on the order parame- 
ter apart from the conventional thermodynamic variables 
and only after minimization with respect to the order 
parameter and substitution of its equilibrium value into 
the original formula one obtains the equilibrium thermo- 
dynamic potential. The order parameter describes the 
structural changes in the material leading to the break- 
ing of symmetry, which is the central idea of this the- 
ory. In the case, of e.g. cubic-tetragonal phase tran- 
sition in BaTiOs this structural change consists mainly 
in shift of the Ti ion with respect to the Ba sublattice. 
Of course, the O ions also displace but to change the 
symmetry from cubic to tetragonal the shift of Ti ions 
is sufficient. Equally, a shift of only O ions without a 
shift of Ti ions would lead to the same symmetry break- 
ing. This illustrates the statement that in the Landau 
theory the physical meaning of the order parameter does 
not matter but only its symmetry properties. For ferro- 
electrics, this is emphasized in the book by Strukov and 
Levanyuk 61 and is correct while the effects of the depo- 
larizing field are not considered. Quite often, it makes 
no difference whether to use P or D within the region 



of applicability of the "orthodox" Landau theory, i.e. for 
small values of the order parameter and not far from the 
phase transition. It is known, in particular, that in ferro- 
electrics normally P ^ E, i.eAirP ~ D. However, even 
close to the phase transition one has to take into account 
that the physical order parameter represents structural 
changes with respect to which the system loses its sta- 
bility at the transition, while considering effects of the 
depolarizing field. In the case of proper ferroelectrics, 
these structural changes can be related to a polarization 
but not necessarily to the full polarization of the crystal. 
At the same time, it is the latter which enters the elec- 
trostatics equations. This can make it necessary to take 
into account the " background dielectric constant" , which 
is an independent phenomenological constant not enter- 
ing the "orthodox" Landau expansion, to be determined 
from the experimental data or from the microscopic the- 
ory. 

We see that much care is needed while applying the 
LGD theory to discuss different phenomena in ferro- 
electrics. The examples where the proper care was not 
taken are too numerous to be discussed here. One can say 
only that the book by Strukov and Levanyuk^ can be im- 
proved in this regard. As an illustrative example, recall 
the boundary condition for the order parameter imposed 
by Ivanchik- 2 who considered a homogeneously polar- 
ized ferroelectric slab with a depolarizing field screened 
by the charge carriers of the material. The electrostat- 
ics dictates that the value of D n , should be zero at the 
surface. Since Ivanchik considered D n as the order pa- 
rameter entering the LGD functional he was forced to 
conclude that the order parameter is zero at the surface. 
But what does it mean? That the Ti or O ions in the 
surface layer are forced to stay in the same positions as 
in the cubic phase? Nonsense. These ions can be shifted, 
of course, but in such a way that at the surface D n = 0. 
This is clearly seen if one considers P n as the order pa- 
rameter. One concludes that at the surface E n = — 47rP n 
but these quantities are not necessary zero and what are 
their values is a separate question. If one takes into ac- 
count the background dielectric constant (e^) one obtains 
E n = — 47rP„/efc, which does not change the qualitative 
result but may be important for numerical estimations. 
Batra and Silverman^, (see Historical note) were right in 
their criticism of Ivanchik's boundary condition in spite 
of considering a case where this condition is not necessary 
and where using either P n or D n as the order parameter 
makes no difference, and it is also not necessary to take 
into account e^. 

There are several equivalent forms of writing down 
the free energy of ferroelectric taking into account elec- 
tric field created by the polarization itself and by exter- 
nal voltage source in the cases where there is no other 
source of the field rather than the ferroelectric polariza- 
tion (like in a non-electroded ferroelectric slab) or there 
is no electric field outside of the ferroelectric and the 
voltage source (like in a ferroelectric slab with the ideal 
metallic electrodes). We shall not discuss these forms 
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here because we are interested in a more general case. 
Let us mention only that a careless generalization to this 
case of one of the above-mentioned formulas leads quite 
often to errors, which could be avoided if the authors 
were to compare at least some of their results with those 
obtained from the system of equations of state and equa- 
tions of electrostatics. Obviously, the equations of state 
usually generate much less controversy than the forms of 
the free energy. 

In the present review, we use nearly the same formula 
for the free energy as Chensky and TarasenkoiA, with 
the only difference in writing a noncritical polarization 
term. They did not discuss the origin of their formula 
but we shall explain why we think it is quite natural. 
Conceptually, it follows the same line as, e.g. the cal- 
culation of energy of collective excitations accompanied 
by long-range (electric or magnetic) fields^. First, one 
considers these excitations without any long-range fields 
as if these fields were compensated by external charges 
or currents, the so-called "mechanical excitons". Then, 
one adds the energy of the long-range field created by the 
exciton without an account for the changes, which such 
fields could produce in the exciton itself, i.e. as if the 
fields were created in vacuum by the bound charges or 
the Ampere currents of the " mechanical exciton" . In our 
case, the "mechanical" part for the ferroelectric energy is 
given by Eq. |T]) , which does not contain any electric field 
and P in this formula has to be considered as a purely 
structural parameter. There is of course also a "mechan- 
ical" noncritical part of the dead layer. The both have 
the form of the first term of the LGD formula (with an- 
other coefficient, of course), because we do not take into 
account neither nonlinearities (the second term) nor the 
spatial dispersion (the third and the fourth terms) in the 
noncritical dielectric properties and in those of the dead 
layer. It is shown in the text, see Eq. (|68|) . that the sum of 
the two energies, the "mechanical" and the "purely elec- 
tric" ones, gives what is used to be called the "energy of 
electric field in a dielectric" . If we had considered explic- 
itly a real metallic electrode, which we avoided for sim- 
plicity sake, we would be obliged to take into account not 
only the energy associated with the displacements of its 
atoms and deformations of the filled electronic shells (the 
"background dielectric constant"), but also an energy as- 
sociated with redistribution of electrons in the conduc- 
tion band, while calculating the " mechanical" energy of 
the electrode. This part is quite often overlooked but a 
detailed discussion of this question would take too much 
space to do it here. 



APPENDIX D: LANDAU PARAMETERS FOR 
STRAINED BATIO3 

The Landau parameters entering the equation of state 
of perovskite BaTi03 film on SrRuOa/SrTiOa entering 
Eq.® arc found from the parameters for the bulk BTO 24 



under conditions of homogeneous in-plane misfit strain 
u m = (a - a )/a = (6 - b Q )/b Q = -0.022: 



A = 2al B = 4a* 33 , 



C = 60111, F = 8oun, (Dl) 



where 



ai 



2Qi 



Sll + «12 ' 



^33 



an + 



Q12 



Sll + S12 



(D2) 



The dielectric function in direction along the 
film ej_ = 1 + (2eoa*) , where a* = a\ — 
u m (Qu + Q12) I [s\\ + Si2)^. When the strain drives 
transition second order from the weak first order, one 
can approximate the dielectric function in a direction 
along the easy axis z and P z = P3 (perpendicular to 
the film plane) as e c = 1 + (4e lagl) -1 . The bulk 
BTO parameters appearing above are (in SI units): 
ai = 4.124 x 10 5 (T - 115) CT 2 m 2 N, where T is the 
temperature in °C, an = -2.097x 10 s C~ 4 m 6 N, a m = 
1.294 x 10 9 C- 6 m 10 N, a mi = 3.863 x 10 10 C- 8 m 14 N, 
Qn = 0.10, Qi 2 = -0.034, Q 44 = 0.029 (C^m 2 ) 28 . 
The components of the compliance tensor for BTO 



are sn = 8.3 x 10" 



S12 



= -2.7 x 10~ 12 , and 



S44 = 9.24 x 10~ 12 (in units of m 2 N _1 ). 

Here in the text we are using the standard CGS units 
instead of the SI units, and the quantity we need to make 
estimates is the (dimensionless in CGS) first Landau co- 
efficient A = eoAgi = —0.0054, e± = 218 at room tem- 
perature. The coefficient Agi is found from (|D1[> . (|D2|) 
and other relations given above. 



APPENDIX E: GRADIENT TERM FOR BATIO3 

The gradient coefficient D for BaTi03 was estimated 
from the neutron scattering data 2 ^ in the following way. 
The experimental data gives the dependence of the soft 
mode frequency on the wave vector for the directions 
perpendicular to (100). In terms of the present paper, 
the soft mode frequency in the paraelectric phase can be 
obtained from the polarization dynamic equation: 



d 2 P k 

1 



(A + Dk 2 ) P k = 0, 



(El) 



where m is a coefficient which will be determined, and 
k is the wave vector in the plane perpendicular to (100). 
Then, 



to 2 (jfc) 



A + Dk 2 



(E2) 



We have found the coefficient to from the value for oj (0) , 
and then estimated from the rest of the curve V^o-Dsi = 
0.17A (Dsi is the gradient term in SI units, the same 
estimate as we gave in ReL 1 ^). In the CGS units, used 
in the present paper, the value is -J D/4tt = 0.17A. 
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